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ABSTRACT 

A full particle simulation study is carried out on the electron acceleration at a collisionless, rela- 
tively low Alfven Mach number (M^ — 5), perpendicular shock. Recent self-consistent hybrid shock 
simulations have demonstrated that the shock front of perpendicular shocks has a dynamic rippled 
character along the shock surface of low-Mach-number perpendicular shocks. In this paper, the effect 
of the rippling of perpendicular shocks on the electron acceleration is examined by means of large-scale 
(ion-scale) two-dimensional full particle simulations. It has been shown that a large-amplitude electric 
field is excited at the shock front in association with the ion-scale rippling, and that reflected ions are 
accelerated upstream at a localized region where the shock-normal electric field of the rippled structure 
is polarized upstream. The current-driven instability caused by the highly-accelerated reflected ions 
has a high growth rate to large-amplitude electrostatic waves. Energetic electrons are then generated 
by the large-amplitude electrostatic waves via electron surfing acceleration at the leading edge of the 
shock transition region. The present result suggests that the electron surfing acceleration is also a 
common feature at low-Mach-number perpendicular collisionless shocks. 
Subject headings: acceleration of particles — plasmas — shock waves — 



1. INTRODUCTION 

Collisionless shocks are universal processes in space 
and are observed in laboratory, astrophysical, and space 
plasmas, including astrophysical jets, an interstellar 
medium, the heliosphere, and planetary magnetospheres. 
Particle acceleration is a common feature of at collision- 
less shocks, but is still a major unresolved issue. The 
most plausible mechanism of the acceleration is the dif- 
fusive shock acceleration (DSA), which explains broad- 
band power-law spectrum with a n index around 2 (Drur^ 
fl98l iBlandford &: Eichlerj ll9871. Before the DSA phase 
in which electrons cross the shock front many times, they 
have to be pre-accelerated by unknown injection mecha- 
nism. Such "injection problem" is still unresolved. 

One of the possible injection mechanisms for electrons 
is the electron surfing acceleration by electrostatic fields 
(iKatsouleas fc D awsoni 119831: IShimada & Hoshinol l2000|; 
iHoshino fc Shim ada 2002). Shimada & Hoshino (2000) 
performed a one-dimensional (ID) electromagnetic full 
particle simulation and found the formation of large- 
amplitude electrostatic solitary structures during the 
cyclic reformation of a high-Mach-number perpendicu- 
lar shock in a low-beta and weakly magnetized plasma. 
The coherent solitary structures trap electrons in their 
electrostatic potential well, resulting in significant surf- 
ing acceleration of electrons. 

So far, the electron surfing acceleration has been 
argued for high-Mach-num ber shocks, and stud- 
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ied by many auth or s (e.g ., . 

[2003 : [Schmitz et all 12002'allb 
2007f l. A typical astrophysical 



example of such high-Mach-number shocks is the young 
supernova remnant (SNR), which is thought to be an 
accelerator of high-energy electrons, emitting bright syn- 



chrotron radiation in radio and X-ray bands at shocks 
with Mach number of M a ^, 1 0^ (e.g., Kovama et al.l 
119951: iBamba et al] l2003l boOSD . However, in many 
astrophysical environments, the electron acceleration at 
lower-Mach-number {Ma ^ 10^) shocks is inferred from 
observations and/or expected theoretically. Possible ex- 
amples are old SN Rs or old SNRs interacting with giant 
molecular clouds ( jClaussen et al. 1997; Chevalier 1991; 
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Bvkov et all 120001 : [Yamazaki et al. 2006.: 

20061) ■ merging galaxy clusters (jMarkevitch et al.li^v^v^^i. 
AGN outbursts at the galaxv clusters (jFuiita et al.1 
[20f)l . and so on. The electron acceleration is also 
observed at Earth's b ow sh ock even with a low Mach 
number of Ma ^ 6.4 (jOka et al. ,200 6). Hence it IS in- 
teresting to ask whether the electron surfing acceleration 
can occur at lower-Mach-number shocks or not. 

It has been demonstrated by one-dimensional full par- 
ticle simulations that the electron surfing acceleration 
takes place at a relatively-high-Mach-number shock of 
Ma ^10 but not at a very-low-Mach-number shock of 
Ma - 3 fe.g.. IShimada fc HoshiiMll2000t [Schmitz et all 
l2002a| ). ^ However, these previous works on the electron 
surfing acceleration are mainly based on ID simulations, 
in which configuration of shock magnetic fields cannot 
be modified. 

It is well known from resuts of two-dimensional hybrid- 
code simulations that there exist large-amplitude fluc- 

^ Note that in the full particle simulations, a reduced ion-to- 
electron mass ratio of mi/nie = 20-25, is usually adopted for com- 
putational efficiency. In general, different mass ratio sometimes 
changes the simulation results qualitatively, and the simulation re- 
sults with a reduced mass ratio cannot be directly compared with 
astrophysical phenomena or results of hybrid simulations. There- 
fore, to avoid a confusion, we use a notation Ma for the Mach 
number of shocks with a reduced mass ratio, while we use another 
notation A4a for the real mass ratio. 
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tuations in the magnetic field and density of the shock 
transitio n region of lower-Mach -number shocks {A4a = 
5 ~ 10) (jWinske fc Questlll988D . These fluctuations at 
the shock surfac e have been analyzed in terms of tur - 
bulent " r ipples" (|Lowe fc Bu rgess '20031 lBurgessl[2006ah . 
iBurgessI ()2006b[ ) has demonstrated, using a combina- 
tion of self-consistent hybrid simulation and test par- 
ticle calculation, that energetic electrons are produced 
by both magnetic mirroring by the rippled structure 
and by stochastic acceleration by magnetic fluctuations 
keeps particles within the shock transition region. How- 
ever, electron-scale microscopic instabilities, i.e., current 
driven instabilities are neglected in the combination of 
hybrid simulation and test-particle calculation. Hence 
multi-dimensional full particle simulations are necessary 
for studying the effect of ion-scale rippling of a perpen- 
dicular shock on the electron surfing acceleration. 

The purpose of this paper is to examine the effect of 
the dynamic rippled structuring of the shock front on 
electron acceleration processes by the first-principle full 
particle simulation. In order to take into account the 
rippling of a perpendicular shock, simulation domain is 
taken to be larger than the ion inertial length. 

The paper is organized as follows. Section 2 describes 
the model and the parameters of the full particle simula- 
tion. Section 3 demonstrates the ion-scale structure of a 
perpendicular shock found in the full particle simulations 
and the resulting electron acceleration. Finally, section 
4 gives summary of this paper. 

2. FULL PARTICLE SIMULATION 

In this paper, a coUisionles s shock is excited with 
the "relaxation method" fe.g.. iLerov et al.lll98lL 119821 : 
lUmeda fc Yamazakil l200l lUmeda~et al.ll2008| ) in which 
the simulation domain is taken in the rest frame of the 
excited shock. Generally speaking, it is not easy to per- 
form a large-scale (ion-scale) multi-dimensional full par- 
ticle simulations of coUisionless shocks even with present- 
day supercomputers. This is because a shock wave ex- 
cited by conventional method becomes unsteady relative 
to the simulation domain. Very recently, however, a 
two-dimensional shock-rest-frame model has been suc- 
cessfully developed bv. Umeda et al.i (|2008.1 . which is im- 
portant to be able to follow the shock for a long time 
with a limited computer resource. 

We use a 2D electromagnetic particle code (jUmedal 
|2004[ ). in which the full set of Maxwell's equations and 
the relativistic equation of motion for individual elec- 
trons and ions are solved in a self-consistent mannar. 
The continuity equation for charge is also solved to com- 
pute the exact c urrent density given by the motion of 
charged particles (jUmeda et al.ll2003[ ). 

The initial state consists of two uniform regions sepa- 
rated by a discontinuity. In the upstream region that is 
taken in the left hand side of the simulation domain, 
electrons and ions are distributed uniformly in space 
and are given random velocities {vx,Vy,Vz) to approx- 
imate shifted Maxwellian momentum distributions with 
the drift velocity u^i, number density ni = eQ-nieUj^^il , 

isotropic temperatures T^i = nieV^^i and Tn = rriiV^j^i, 
where m, e, ujp and Vt are the mass, charge, plasma fre- 
quency and thermal velocity, respectively. Subscripts "1" 
and "2" denote "upstream" and "downstream" , respec- 
tively. The upstream magnetic field ByQi = —m^ujcei/e 



is also assumed to be uniform, where Wc is the cyclotron 
frequency (with sign included). The downstream region 
taken in the right-hand side of the simulation domain is 
prepared similarly with the drift velocity 7/3,2, density 712, 
isotropic temperatures Tf,2 and Ti2, and magnetic field 

We take the simulation domain in the x-y plane and 
assume a perpendicular shock (i.e., B^q — 0). Since the 
ambient magnetic field is taken in the y direction, free 
motion of particles along the ambient magnetic field is 
taken into account. As a motional electric field, a uni- 
form external electric field Ezo = —u^iByoi = —Ux2By02 
is applied in both upstream and downstream regions, so 
that both electrons and ions drift in the x direction. At 
the left boundary of the simulation domain in the x direc- 
tion, we inject plasmas with the same quantities as those 
in the upstream region, while plasmas with the same 
quantities as those in the downstream region are also 
injected from the right boundary in the x direction. We 
adopted absorbing boundaries to suppress non-physical 
reflection of electromagnetic waves at both ends of sim - 
ulation domain in the x direction (jUmeda et al.|[200ll ). 
while the periodic boundaries are imposed in the y di- 
rection. 

The initial downstream quantities are given by solv- 
ing the shock jump conditions for a magnetized two-fluid 
isotropic plas ma consisting of electrons and ions (e.g., 
lHudsonl[l970h . 

^lel'^xl = Ujl^2Uxk\) 
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where Ts = nisvjg. In order to determine a unique ini- 
tial downstream state, we need given upstream quantities 
Uxi, ojpei, LJcel, Vtel, and Tn/Tpi and an additional pa- 
rameter. We assume a low-beta and weakly-magnetized 
plasma such that /?ei — Pn = 0.125 and LUeei/(^pei — 
—0.1 in the upstream region. We also use a reduced ion- 
to-electron mass ratio rrii/me = 25 for computational 
efficiency. The light speed c/vtei = 40.0 and the bulk 
flow velocity of the upstream plasma Uxi/vtei — 4.0 
are also assumed. Then, the Alfven Mach number is 
calculated as Ma = (uxi/c) \ujpei / tOcei | / irie = 5.0. 
The ion-to-electron temperature ratio in the upstream 
region is given as Tn/Tei = 1.0. In this study, down- 
stream ion-to-electron temperature ratio Ti2/Te2 = 8.0 
is also assumed as another initial parameter to obtain 
the unique downstream quantities by solving the above 
four equations, iOpe2/^pei — 1.8372, U!ce2/^pei = 0.3375, 
Ux2/vtei = 1.1851, and Vte2/vtei = 2.6393. 

In this study, we perform two runs with different sizes 
of the simulation domain. We use NxXNy — 2048 x 1024 
cells for the upstream region and Nx x Ny — 2048 x 1024 
cells for the downstream region, respectively, in Run A. 
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The grid spacing and time step of the present simulation 
are /^.x/Xoei — 1-0 and LOpei/Si — 0.0125, respectively. 
Here Xoei is the electron Debye length upstream. Thus 
the total size of the simulation domain is 10.24Ai x 5.12Ai 
which is long enough to include the ion-scale rippled 
structure, where Xi = c/ujpii{— 2QQXDei) is the ion iner- 
tial length. In Run B, we use iVj, x iV^ = 2048 x 128 cells 
for the upstream region and x Ny = 2048 x 128 cells 
for the downstream region, respectively. Thus the the 
total size of the simulation domain is 10.24Ai x 0.64Ai, 
in which ion-scale processes along the ambient magnetic 
field is neglected. We used 16 pairs of electrons and ions 
per cell in the upstream region and 64 pairs of electrons 
and ions per cell in the downstream region, respectively, 
at the initial state. 

3. RESULT 

Figure 1 shows the transverse magnetic field By as a 
function of position x and time t averaged over the y di- 
rection for Run A. The position and time are renormal- 
ized by the ion inertial length A^ and the ion cyclotron 
period l/wcii, respectively. The magnitude is normal- 
ized by the initial upstream magnetic field ByQi. In the 
present shock-rest-frame model, a shock wave is excited 
by the relaxation of the two plasmas with different quan- 
tities. Figure 1 shows that the shock front appears and 
disappears at a timescale of the downstream ion gyro- 
period, which corresponds to the cyclic reformation of 
perpendicular shocks. Since the initial state is given by 
the shock jump conditions for a two- fluid plasma consist- 
ing of electrons and ions, the excited shock is "almost" 
at rest in the simulation domain. 

We have also performed a ID simulation with the same 
parameter as the 2D simulations and confirmed that the 
period of cyclic reformation in all simulation runs (Runs 
A, B, and ID) is the same. However, the peak ampli- 
tude of the overshooting magnetic field in Run A is about 
6.0i3yoii while this is about 4.5-Byoi in Run B and in the 
ID run. We have also confirmed that the cyclic reforma- 
tion becomes less significant for uJdit > 7 in Run A, while 
the cyclic reformation continues for a long time in Run 
B and in the ID run. The present result is in agreement 
with the recent 2D full-particle simulations of a low-Mach 
number shock, in which it has been demonstrated that 
the perpendicular shock evolves from the cyclic reforma- 
tion phase to the "nonlinear whistler" phase i n 2D simu- 
lations with the in-plane m agnetic field (|HeUinger et al.l 
[20071: iLembege et al.l[20M) . 

In the present study we focus on the particle accelera- 
tion in the cyclic reformation phase. Figure 2 shows the 
rippled structure of the perpendicular shock at uJcut = 
6.256. The top panel shows a gray-scale map of the mag- 
netic field magnitude By, and the bottom panel shows a 
gray-scale map of the electric field magnitude around 
the shock transition region. The magnitude of magnetic 
field is normalized by the initial upstream magnetic field 
ByQi, while the magnitude of electric field is normalized 
by the motional electric field Ezo- 

We found a strong fluctuation in the magnetic field 
component By at the shock surface. The amplitude 



of the fluctuation is estimated as ^ 1.5B 



t/Ol, 



which is 



larger than the magnitude of the upstream magnetic 
field. The wavelength of the fluctuation is several ion 
inertial length. This strong fluctuation has been ana- 
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Fig. 1. — Transverse magnetic field By as a function of posi- 
tion X and time t averaged over the y direction for Run A. The 
position and time are normaHzed by and l/ujcn, respectively. 
The magnitude is normalized by the initial upstream magnetic field 

ByOl- 



lyzed by iLowe &: Burgess! (|2003l ) in terms of ripples at 
the shock overshoot. 

We computed a numerical frequency-wavenumber 
spectrum by taking Fourier transformation of the B^ 
magnetic field component (not shown) and obtained a 
similar spectrum to the result bv lLowe" fc Burgess! (|2003D . 
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Fig. 2. — The rippled structure around the transition region 
of the perpendicular shock at utciit = 6.256 for Run A. The mag- 
netic field magnitude By (top) and the electric field magnitude Ex 
(bottom). The magnitude of magnetic field and the electric field 
is normalized by the initial upstream magnetic field Byoi i and the 
motional electric field Ezo, respectively. 

The generation mechanism of the ripples is thought to 
be the ion perpendicular temperature anisotropy in the 
shock transition region, which drives Alfven ion cyclotron 
(L-mode) or mirror mode waves. Note that there is not 
any rippled structure in Run B because the size of the 
simulation domain in the magnetic field (y) direction is 
shorter than the ion inertial length. The overshoot mag- 
netic field By is uniform along the magnetic field in Run 
B and its magnitude is about 4.5i3yoii which is equal to 
the average overshoot magnetic field in Run A. 

Around the shock front {x/Xi = —0.5 ~ 0), we also 
found a strong negative electric field with a magnitude 
of ~ — 4.0£'zo (see the bottom panel of Figure 2). This 
strong negative electric field at the shock front refiects 
ions upstream. Such a negative electric field is also found 
in Run B, but its magnitude is much smaller, about 
Ex/ ~ —2.0. There is a good correlation between 
the structure of By and E^ at the shock front as seen in 
Figure 2, implying that the strong negative electric field 
is associated with the rippled structure. 

The effect of the acceleration of refiected ions is ana- 
lyzed in Figure 3. The electric field magnitude at 
ijJciit = 6.575 for Run A is shown in Figure 3a. At 



[x/Xi.y/Xi) ~ (-1.0,1.5) and (-1.0,3.5), we found wave 
structures with a short wavelength, which corresponds 
to electrostatic waves excited by the current-driven in- 
stability due to reflected ions. An interesting result here 
is that the electrostatic waves are not excited uniformly 
in the y direction but are excited in localized regions. It 
is suggested that the localized excitation of the electro- 
static waves is due to the strong negative shock-normal 
electric field at the shock front that accelerates refiected 
ions upstream. 

The strong negative shock-normal electric field is as- 
sociated with the rippled strucuture. Figures 3c and 
3d show the mangetic field components By and Bz at 
Wciii — 6.575 for Run A. We have analyzed the Hall elec- 
tric field (J X B)/en and found that the term of JzBy/en 
is dominant. In other words, the strong negative shock- 
normal electric field is due to the magnetic pressure gra- 
dient force of the By componet around x/Xi ~ —0.5 aris- 
ing from the rippled structure. Here, the excitation of 
electrostatic waves is discussed in terms of the particle 
distributions shown in Figures 3b and 3e. 

Figure 3b shows the ion density rii normalized by the 
upstream density ni. Around y/Xi ~ 0.2, 2.0, and 3.2, 
we found the enhancement of the ion density that exceeds 
15ni. In these regions, the magnitude of the magnetic 
field becomes lower (see By in Figure 3c) to keep the 
pressure balance. 

The left panels of Figure 3e show y — Vx phase-space 
plots of ions and electrons. Around y/Xi ^ 0, 1.5, and 
3.5, we found ion components with negative velocity, 
which correspond to reflected ions. The electrostatic 
waves around y/Xi 1.5 and 3.5 are excited by the local- 
ized reflected ions. The patchy reflected ion beams are 
formed by the localized negative shock-normal electric 
fleld around x/Xi = —0.5. 

Around y/Xi ~ 0, despite a strong shock-normal elec- 
tric fleld, electrostatic waves are not strongly excited, 
either. This is because of small number of ions. We 
need a strong ion current in the —x direction to enhance 
the current-driven instability. That is, the velocity of 
reflected ions should be high and/or the number of re- 
flected ions should be high to excite the patchy electro- 
static waves. 

In the y — Vz phase-space plots in the right panels, non- 
thermal components of electrons are found at y/Xi ~ 1.5, 
where electrostatic waves are strongly enhanced. The ex- 
istence of nonthermal components in Vz is the evidence 
of the electron surflng acceleration at low-Mach-number 
perpendicular shock with Ma = 5. The electron surf- 
ing acceleration takes place in a localized region where 
a negative shock-normal current due to reflected ions is 
strongly enhanced. 

Figure 4 shows energy distribution functions of elec- 
trons in the downstream region for Runs A and B. The 
solid line shows the result of Run A, and the dashed 
line shows the result of Run B. Firstly, we found non- 
thermal electrons in Run A while there is no nonthermal 
electrons in Run B. Secondly, the downstream temper- 
ature in Runs A and B is very different, indicating the 
existence of additional electron energization mechanism 
in Run A, i.e., rippling of the perpendicular shock dis- 
cussed by [Burgess! (|2006bi) . We have also performed a 
ID simulation run and conflrmed that the electron en- 
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Fig. 3. — (a) The electric field magnitude Ex at uidit = 6.575 for Run A. (b) The corresponding ion density n^. (c) The corresponding 
magnetic component By . (d) The corresponding magnetic component . (e) The corresponding reduced phase-space distribution functions 
in y — Vx and y — Vz, at x/\i = —1.0. The density is normalized by the upstream density ni. The velocity is normalized by the upstream 
bulk velocity Uxi- 



6 



Umeda, Yamao & Yamazaki 



Electron Energy Spectra 



XI 

E 
=1 
-z. 

c 

o 

I 

LU 



10' 



10' 



10' 



10' 



10 °0 





■ Run A 


V 


RunB ■ 










* ■■ ^ 
1 

\ 

> * 





10 



20 
V„2 / U, 



30 



40 



Fig. 4. — Energy distribution functions of electrons in tlie down- 
stream region. Tlic solid line shows the result of Run A and the 
dashed line shows the result of Run B. The energy is normalized by 
the upstream bulk energy of electrons. The dotted line indicates a 
Maxwellian distribution for the downstream of Run A. 

ergy spectrum is essentially the same as that in Run B. 

The electron energy distribution for Run B (dashed 
line) corresponds to the Maxwellian distribution func- 
tion with thermal velocity oivte ~ 2.3, which can be ex- 
plained by the plasma heating by magnetic compression. 
On the other hand, the thermal component of electron 
energy distribution for Run A (dotted line) corresponds 
to the Maxwellian distribution function with thermal ve- 
locity of Vte ~ 3.1. In addition, the maximum energy of 
electrons reaches ^ 30u^]^ in Run A. 

The rippled structure can further energize electrons 
as discussed by lBurge"ssl ()200 6b). which can explain the 
Maxwellian distribution for Run A. In the shock foot re- 
gion (xjXi = —1 ~ 2), the magnitude of the magnetic 
field is By ^ 1.5Byoi. The average amplitude of electro- 
static waves excited by the current-driven instability is 
Ex ~ 'iEzo- Thus the maximum energy of non-thermal 
electrons becomes about iu^i ~< + vl >= \Ex/By\^ 
via the electron surfing acceleration as indicated from the 
y — Vx and y ~ Vz phase-space plots shown in Figure 3e. 
Thus there is three effects for energization of non-thermal 
electrons: The first is the electron surfing acceleration; 
The second is the magnetic compression at the shock 
overshoot (e.g..[Shimada fc Hoshinoll200Cll : [Schmitz et al.l 
l2002aHU^da fc Yamazakill2006[ ). The third is the scat- 
tering by ion-scale rippled structure (|Burgessll2006b[ ) . To 
reach the maximum energy of SOu^i all the three ef- 
fects are needed. 

Both in Run B and in the ID run, the current- 
driven instability due to refiected ions excites electro- 
static waves with amplitude of ~ 2.0Ezo. However, 
this amplitude is not enough for the electron surf- 
ing acceleration. By contrast, the amplitude of the 
locally-excited electrostatic wave exceeds ~ i.OEzo in 
Run A as shown in Figure 3a, and energetic elec- 
trons are generated by the electron surfing accelera- 
tion. This result suggests that microscopic current- 
driven instability enhanced by the ion-scale rippled struc- 
ture of the low-Mach-number perpendicular shock gen- 
erates nonthermal electrons, as observed in moderate- 
Mach-number perpendicular shocks of Ma = 10 — 
20 (e.g. Shimada fc Hoshino 2000: Hos hino fc Shimadal 
I2OO2I: iSchmitz et al.H2002aibl : lAmano fc Hoshinoll2007ir ^ 



4. SUMMARY AND DISCUSSION 

We have studied electron acceleration at a low-Mach- 
number perpendicular coUisionless shock by perform- 
ing two-dimensional full particle simulations. In or- 
der to take into accou nt the effect of the rippling of 
a perpe ndicular shock ijLowe fc Burgess! 120031 : iBurgessI 
l2006allbl ). the simulation domain is taken to be larger 
than th e ion inertial scale by using the shock-rest-frame 
model (jUmeda et al.ll2008D . In the previous works, the 
electron acceleration by electron-scale microscopic insta- 
bilities at low-Mach-number perpendicular shocks has 
not focused on because the e lectron surfing acceleration 
by elec trostati c wav es _(e.g., Shimad a fc Hoshinol 120001 : 
iHoshino fc Shimadal l200§ ) . which is one of the possible 
injection mechanisms for electrons, is thought to be ef- 
fective only in high-Mach-number perpendicular shocks. 
The present result suggests, by contrast, that the elec- 
tron surfing acceleration is also a common feature at a 
perpendicular shock with — 5. 

The mechanism for generation of high-energy electrons 
is quite simple, and is summarized as follows: The per- 
pendicular shock forms the rippled structures by ion 
temperature anisotropy; The rippled structures excite a 
strong electric field component in the shock-normal di- 
rection; Reflected ions in the shock transition region is 
strongly accelerated upstream by the electric field com- 
ponent of rippled structures; The strongly-accelerated re- 
flected ions give a high growth rate of a current-driven 
instability to large-amplitude electrostatic waves in a 
localized region; Energetic electrons are generated by 
the electrostatic waves via surfing acceleration. Finally, 
high-energy non-thermal electrons are generated by both 
magnetic compression and scattering by ion-scale rippled 
structure in the shock transition region. 

It has been demonstrated that multi-dimensionality 
sometimes weakens the electron surfing acceleration in 
the 2D simulations of uniform plasma models (e.g., 
iDieckman n fc Shuklal l2006l:IOhira fc Takahara"2007|iaiid 
self-c o nsistent shock mo dels (e.g., Dicckman n et al.l 
2008t lUmeda et al.| I2008D . Very recently, however. 



Amano fc Hoshinol (j2008[ ) have found the electron surfing 



acceleration in 2D simulation of a perpendicular shock 
with out-of-plane magnetic fleld. In their 2D simulation, 
the simulation domain is taken to be much larger than 
the ion inertial length along the shock surface, and ki- 
netic effects of ions are fully included. The present result 
is another example showing the effect of ion kinetics to 
the electron surflng acceleration. It has been conflrmed 
that the cross-scale coupling between an ion-scale meso- 
scopic instability and an electron-scale microscopic in- 
stability is important. Hence, a large-scale full particle 
simulation would be essential for studies of the electron 
acceleration at coUisionless shocks. 

Finally let us discuss on the reduced mass ratio. In 
the present simulation with Ma — 5, LOpei/ucei = 10, 
and /3 = 0.125, the Buneman-type mode becomes un- 
stable at LUuHR — ^pei when we use a reduced mass 
ratio of mi/rUe — 25. When we use the real mass ra- 
tio of mi/rrie — 1836, the thermal velocity of upstream 
electrons becomes about 8.6 times as large as the case 
of mi/me = 25, the self-reformation process is sup- 
pressed and the Buneman-t ype mode is also stabihzed 
(jScholer fc Matsukivol l2004l ). To make the Buneman- 
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type mode unstable the electron thermal velocity must 
be smaller than the relative velocity between incoming 
electrons and reflected ions. Thus the present simulation 
result can be applied to a much lower beta of /? < 0.03. 

The authors are grateful to Y. Ohira and Y. Fujita for 
helpful comments. The computer simulations were car- 
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rative computational research project at STEL, Nagoya 
University, and YITP, Kyoto University. This work was 
supported by Fujii research grant from Hiroshima Uni- 
versity, Grant-in-aid 18740153, 19047004 (R. Y.), and in 
part 17GS0208 (T. U.) from MEXT of Japan. 



REFERENCES 



Amano, T., & Hoshino, M. 2007, ApJ, 661, 190 
Amano, T., & Hoshino, M. 2008, ApJ, in press 
Bamba, A., Yamazaki, R., Ueno, M., and Koyama, K. 2003, ApJ, 
589, 827 

Bamba, A., et al. 2005, ApJ, 621, 793 

Blandford, R. D. & Eichler, D. 1987, Phys. Rep. 154, 1 

Brogan, C. L. et al. 2006, ApJ, 639, L25 

Burgess, D. 2006a, J. Geophys. Res., Ill, A10210 

Burgess, D. 2006b, ApJ, 653, 316 

Bykov, A. M. et al. 2000, ApJ, 538, 203 

Claussen, M. J. et al. 1997, ApJ, 489, 143 

Chevalier, R. A. 1999, ApJ, 511, 798 

Dieckmann, M., & Shukla, P. K. 2006, Plasma Phys. ControU 

Fusion, 48, 1515. 
Dieckmann, M. E., MeU, A., Shukla, P. K., Drury, L. O. C. & 

Mastichiadis, A. 2008, Plasma Phys. Controll Fusion, 50, 065020. 
Drury L. O'C. 1983, Rep. Prog. Phys. 46, 973 
Fujita, Y. et al. 2007, ApJ, 663, L61 

Hellinger, P., Travnicek, P., Lembege, B., & Savoini P. 2007, 

Geophys. Res. Lett., 34, L14109 
Hoshino, M. & Shimada, N. 2002, ApJ, 572, 880 
Hudson, P. D. 1970, Planet. Space Sci., 18, 1611 
Katsouleas, T., & Dawson, J. M. 1983, Phys. Rev. Lett., 51, 392 
Koyama, K. et al. 1995, Nature, 378, 255 

Lee, R. E., Chapman, S. C, & Dendy, R. O. 2004, ApJ, 604, 187 
Lembege, B., Savoini, B., Hellinger, P., & Travnicek, P. M. 2008, 

J. Geophys. Res., in press 
Leroy, M. M., Goodrich, C. C, Winske, D., Wu, C. S. & 

Papadopoulos, K. 1981, Geophys. Res. Lett., 8, 1269 



Leroy, M. M., Winske, D., Goodrich, C. C, Wu, C. S. & 

Papadopoulos, K. 1982, J. Geophys. Res., 87, 5081 
Lowe, R. E., & Burgess, D. 2003, Ann. Geophys. 21, 671 
Markevitch, M. et al. 2002, ApJ, 567, L27 

McClements, K. G., Dieckmann, M. E., Ynnerman, A., Chapman, 

S. C, & Dendy, R. O. 2001, Phys. Rev. Lett., 87, 255002 
Ohira, Y., & Takahara, F. 2007, ApJ, 661, L171 
Oka, M. et al. 2006, Geophys. Res. Lett., 33, L24104 
Schmitz, H., Chapman, S. C, & Dendy, R. O. 2002a, ApJ, 570, 
637 

Schmitz, H., Chapman, S. C, & Dendy, R. O. 2002b, ApJ, 579, 
327 

Scholer, M., & Matsukiyo, S. 2004, Ann. Geophys. 22, 2345 

Shimada N., & Hoshino, M. 2000, ApJ, 543, L67 

Umeda, T. 2004, Ph.D. Thesis, Kyoto University 

Umeda, T., Omura, Y., & Matsumoto, H. 2001, Comput. Phys. 

Commun., 137, 286 
Umeda, T., Omura, Y., Tominaga, T., & Matsumoto, H. 2003, 

Comput. Phys. Commun., 156, 73 
Umeda, T., & Yamazaki, R. 2006, Earth Planets Space, 58, e41 

( arXiviphysics /0607220 ) 
Umeda, T., Yamao, M., & Yamazaki, R. 2008, ApJ, 681, L85 
Winske, D., & Quest, K. B. 1988, J. Geophys. Res., 93, 9681 
Yamazaki, R. et al. 2006, MNRAS, 371, 1975 



